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We study the effects of Coulomb interaction and inter-grain quantum tunnelling in an array 
of metallic grains using the phase-functional approach for temperatures T well below the charging 
energy Ec of individual grains yet large compared to the level spacing in the grains. When the inter- 
grain tunnelling conductance <? 3> 1, the conductivity cr in d dimensions decreases logarithmically 
with temperature {a/ao ~ 1 — ln{gEc/T))~'~'^, while for g ^ 0, the conductivity shows simple 
activated behaviour (cr ~ exp{—Ec/T)). We show, for bare tunnelling conductance <? > 1, that the 
parameter 7 = g{l'-2/ (gn) ln{gEc/T)) determines the competition between charging and tunnelling 
effects. At low enough temperatures in the regime 1 > 7 S> l/yJPEc, a charge is shared among a 
finite number N = ^ [E^/T) / \n{-K /2'yz) of grains, and we find a soft activation behaviour of the 
conductivity, a ~ exp(— 2y (Ec/T) hi(7r/272:)), where z is the effective coordination number of 
a grain. 

PACS numbers: 

I. INTRODUCTION 

Coulomb effects and electron tunnelling as well as various effects of disorder are major themes of mesoscopic physics. 
For a single metallic grain, all these effects have been explored extensively. For an array of normal metal grains, we 
find a peculiar interplay of electron charging and tunnelling effects manifested by a formation of a multi-grain charge 
excitation. 

The novelty of granular metal systems in comparison with disordered metals or semiconductors with impurities 
arises from the presence of additional energy scales - the grain charging energy and intra-grain energy-level 
spacing S. When the temperature is lowered below level spacing S, electrons in the granular metal can propagate 
(diffuse or hop over many grains) coherently just like in a disordered metal. In the incoherent regime T > max {g6, S), 
only charging effects and nearest-neighbour hopping (in the second order of the tunnelling matrix element) are 
relevant, and it is possible to formulate the problem in terms of longitudinal electromagnetic phase (or voltage) 
fluctuations on the grains^i^. For large inter-grain conductance g ^ 1, the conductivity decreases logarithmically with 
temperature independent of dimensionality^'^'^, reminiscent of many experimentaSi^. The essential physics of this 
result was recognised some time ago as washed out Coulomb blockade for the quantum dot coupled to a conductive 
environmentSiifl. For intermediate conductance g ^ 1 (and sufficiently low temperature), we derive within the same 

approach a soft activation conduction {a ^ e~^'^°^'^) by a charge excitation (later referred to as a "puddle") shared 
between many grains due to incoherent tunnelling. In the itinerant (or large-scale diffusion over many grains) regime 
at low temperatures T < max (gd, S) , the granular metal for large conductance 5 ^ 1 is naturally described by 
Altshuler-Aronov theorjiiiii^. 

We employ the phase functional approach suggested over two decades ago by Ambegaokar, Eckern, and Schon 
(AES)5i. The original model described the tunnelling dynamics of granular superconductors, but nowadays this 
approach is increasingly used to study granular metals as wel&iSii^. The AES action Saes for granular metals 
consists of two contributions, Saes = Sc + St, where 







represents charging of the grains, and 



St = 7rg >^ f' dr dr' air -r') sin' (h^^l-^] (2) 
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represents tunnelling between neighbouring grains, — <j)i — (py The kernel a{r) has the form a(r) — 
T^(Re (sin(7rrT-f ie))^^)^. The fields {0i} are electromagnetic phase fluctuations on the grains related to the respective 
potential fluctuations {Vi} through Vi{r) = dr4>i{r). They satisfy bosonic boundary conditions, (pi{r) = ^^r + (j)iir), 
(j)i{r) = (f>i{r + P), where the winding number ki is an integer, and —00 < (t>i{r) < 00. The tunnelling conductance 
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g is related to the inter-grain hopping amphtude ti,i+a through g = 27r|ti4+aP/^^- Conductivity in the AES model 
is a second order (in hopping amplitude) incoherent tunnelling process between neighbouring grains. The elastic 
tunnelling lifetime t of the electron on the grain is t = h/ (gS). The condition, defining the granularity of the material 
and allowing averaging over fermionic intra-grain states, is that the tunnelling lifetime r is much longer than the 
Thouless diffusion time P/D (where D is intra-grain diffusion coefficient, and / is a size of a grain). Another relevant 
condition is the implicit requirement of energy relaxation in the grains. The characteristic times associated with 
these incoherent dissipative processes should be shorter than tunnelling lifetime, consequently coherent combination 
of wavefunctions over grains cannot be written. Moreover, our diagrammatic analysis shows (see also Ref.^'^) that 
higher order processes (|ti,i+a|*, etc.) can be neglected when T > gS. Such a condition allows us to neglect "dress- 
ing" of the tunnelling vertex ti^+a by further tunnelling lines. Therefore, the phase functional approach for granular 
materials can be justified if g6 <C D/P, and the temperature is sufficiently high, T ^ max{d, gS). We shall restrict 
our analysis to this regime. 



II. ANALYSIS OF MODEL 



The AES action shows important qualitative changes in the relevance of large phase fluctuations as the coupling 
g is varied. Consider the metallic phase 5 ^ 1 in Eq.Q. Expanding sin^(0ij(T) — ^ij(r')) in a power series, we 
observe g~^, thus inter-grain phase fluctuations are Gaussian, and suppressed. The charge on an individual 

grain is not a well defined quantity, rather it is shared by the entire system. As g is progressively reduced, the phase 
fluctuations in Eq.jSJl increase until finally one needs to take into account non-zero winding numbers fci ^ 0. In the 
extreme limit of g — > 0, the AES model describes a system of weakly coupled capacitors. The phase fluctuations are 
large, however the charge on an individual grain is well-defined. Conduction now involves exciting a charge which 
results in an activated temperature dependence of conductivity (cr ~ 5 exp(— /3i?c)). Such considerations lead us to 
examine whether for intermediate coupling between the grains, charge could be shared by a finite number of grains. 
This would be an intermediate situation between the extreme cases discussed above. In the remaining part of this 
letter, we choose a diagonal capacitance matrix in Eq.(^, Cy ~ 2E~'^ij' ^'^ keep our analysis simple. 

We describe now the physical picture for the soft activation phase which emerges from our analysis. Putting a single 
electron on an isolated grain costs i?c, while incoherent tunnelling enables the charge to be shared between two or more 
grains. We show below that the parameter 7 = g{l — 2/ (ng) hi(gPEc)) controls the suppression of winding numbers, 
and determines the degree of charge delocalisation. When 7 > 1, the charge is effectively delocalised over the entire 
system (the charging energy is exponentially suppressed). For 7 1, a unit of charge (electron) is shared among a finite 
number of grains. For simplicity we consider two grains, and compare statistical weights associated with the charge 
localised on any single grain Pi ~ exp(—/3 £"(,), and the charge shared between the two grains, P2 ^ 7exp(— /?£'c/2) 
(/3 = l/T). Observe that the charging energy is halved upon hybridisation^^. Since a charge is shared classically 
(incoherently) between two grains, it is equivalent to equal average voltage on the grains, and thus two capacitors 
connected in parallel. The total capacitance is doubled, and the charging energy is halved. If 7 < exp{—l3Ec/2), the 
charge is unlikely to be shared between the two grains. If on the other hand, 7 > exp{—PEc/2), the electron is more 
likely to live on both the grains. Thus the two-grain hybridisation "puddle" optimising the charging and tunnelling 
energies is formed. The optimum number N^, of hybridised grains sharing a single charge is determined by maximising 
^ j^-'^exp{-(3Ec/N), which gives - ^/pKTHT^, hence a cx gPN, exp{-2y^j3KH^r^). This in 
essence is our main result. 

To calculate the conductivity, we use Kubo's formula'^. 
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-n. Also of interest is the tunnelling density of states Vi^s) into the grain i: 
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where H; = (exp(-i(^i(r) - 0i(O)))), and e„ = ^(n + \). 

At this stage we are in a position to understand qualitatively the logarithmic temperature dependence of the 
conductivity for g ^ 1 (derived in Ref.^-^-'^). Since in this regime phase fiuctuations are small, we set all fci = and 
expand St to quartic order in Denoting the Gaussian part of the resulting AES action as the "free" action, and 
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considering the quartic bit as "interaction" , one finds that the interaction renormahses the tunnelhng conductance 



Qren ('^ 



5(l-((/.ij(r)0y(r'))), 
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thus in d dimensions, greniP) ~ ^1 ~ ln(g/3i?c)^ • One infers from Eq.jSj) a similar temperature dependence for 
the conductivity a as it is proportional to the effective tunnelling conductance gren- This result suggests that when 
temperature T falls below To = Ece'^'^'^^3~^\ a transition or crossover into an insulating phase might be expected. 
Note that the parameter 7 becomes smaller than one at temperature Ti = iJ^e^'^^^"^^/^, which is parametrically 
much larger than Tq. We show that Ti marks the onset of soft activation behaviour. 

Consider the AES model, Eq.QJ and Eq.(j2Il. At low temperatures and for 51 > 1, phase fluctuations on the grains 
could be large and non-Gaussian, so we expand the action about the finite winding number phase changes {2t: / j3)kiT 
and residual fluctuations (br^: 
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where 4ii{T) = (/)i(aj„) exp(ici;„r). Since the bare conductance is large, 5 > 1, an expansion to quadratic order in 
the residual fluctuations is justified. The first two terms of Eq.© arise from finite- winding number (non-Gaussian) 
fluctuations, and directly lead to quantisation of charge. The remaining terms in Eq.JBl arise from perturbation 
about the winding numbers. The competition of single-grain charging and hybridisation at low enough temperatures 
can be seen in the partition function Z2 — ^ D(j)iD(j)2 X^fc^ k2 6xp(— S'[fc, (p]) of a simple two-grain system. Since the 
tunnelling term depends only on the phase difference between the two grains, we make a transformation to average 
phase (j)av = {i'l + </*2)/2, and relative phase = (0i — 02)- Integrating out the relative phase gives a winding number 
dependent determinant. This we normalise against the determinant with no winding numbers: 
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For fci2 ^ gPEc (a large g suppresses large values of A:i2), the determinant in Eq.© can be simplified as 
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where we used Stirling's formula for the factorials. This result implies that relative phase fluctuations enhance the 
tendency for phase slips between neighbouring grains. The effective action for the two grain system takes the form 
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where the effective parameter 



1 = 9 
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is smaller than the renormalised tunnelling gren = 5(1 — ^/ {'^dg)^'^{gPEc)) that represents the conductance at large 
bare coupling g at not too low temperatures. The main interest of this paper is the regime 7^1. lfl^7^1/ (pEc), 
most contribution to single charge excitations in the two-grain partition function comes from low winding number 
difference (|A:i2| ~ 0(1)) between the grains. On the other hand, when 7 < l/{/3Ec), large winding number differences 
of the order of PEc become important. Therefore in this regime, the charge is localised on cither of the grains, with 
large charging energy of order Ec- 

It is important to obtain the charge representation by summing over the winding numbers with the use of the 
Poisson formula 

Z-2 = Y.r i^^^a.^^e'^'^S.^— ^[-'^-1. (11) 

{<?»} 

Integrating out X2 will yield the effective environment around xi. 
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There are two qualitatively distinct contributions in Eq. (|12|l . The first two terms represent isolated charging of grain 
2. If 7 were vanishingly small, this would be the only contribution. The last term represents hybridisation of the two 
grains because of quantum tunnelling; the total charge qi + 92 is shared between the grains, and the charging energy 
Ec is reduced to Ec/2. Finally, integration over xi gives the relative weights of the two processes in the partition 
function as Pi « exp[—q\(3Ec) for isolated charging, and P2 ~ ^ exp(— (gi +(72)^/3£'c/2) for charging of the hybridised 
grains. 

The treatment so far considers residual relative phase fluctuations only to gaussian order. In the Appendix, 
we present the results of path integral Monte Carlo calculations to support our basic idea of charge sharing over 
two grains, even at low temperatures, where non-gaussian fluctuations are important. We also confirm that the 
temperature dependence of 7 agrees well with Ea. (|10|l at not too low temperatures. As the temperature is decreased 
further, non-gaussian fluctuations become important. The numerical calculations show that our physical picture, 
that a competition of charging and tunnelling effects determines whether the charge is shared between the two grains 
(with charging energy i?c/2) or localised on a single grain, still remains valid. More precisely, we find that the 
probability of sharing a charge between two grains (here 7) is an algebraically, and not exponentially, small function 
of temperatureiSi. 

Consider now the case of N connected grains. Formally, it is simple enough to demonstrate annexation of a single 
grain into an A^-site puddle. The proof is by induction. Suppose that an A^-site puddle already exists (with statistical 
weight (X 7^^^). Integration over a string of A' — 1 contiguous neighbours of a grain i similarly gives a puddle of size 
A^ with charging energy E^/N , and a weight Pjv ~ (^)^~^ exp(— (gi -!-••• + qn)'^PEc/N). Such an expansion in 7 
only makes sense if 27/7r < 1. For 27/7r > 1, optimum size of the puddle is divergent. Consider the action of a single 
grain coupled to this puddle: 

(2^)^A^fc^ , {27rn%^, 
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In terms of the centre of mass coordinate 

N(j)N + 4>N+l 



N + 1 



(14) 



and relative coordinate 



the action takes the form 
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Integrating out the relative phase renormalises the bare coupling g in a manner similar to that in Ea. (|10|l . 
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Note that the relevant 7Ar,A'+i determining annexation of a single grain into an N-site puddle is not too different 
from 7 for a two grain system obtained in Ea. (|10|l . Accordingly, the condition for suppression of large winding 
number difference changes from 7 ^ \/{[3Ec) for two grains to the condition 7 ^ N/{PEc) for N grains. Performing 
the summation over fcjv and fcjv+i in Ea. ljlGI) using the Poisson summation formula again yields two terms that 
correspond to separate charging of the puddle and grain, and charging of the larger (N+l)-site puddle. The criterion 
for annexation is 

^ME+l exp[-(3E,/iN + 1)] > exp[-f3EjN]. (18) 

TT 

So far we have obtained the effective environment of a site i by integrating out a sequence of iV — 1 contiguous 
neighbours. Integrating over such 'strings' is somewhat different from the actual requirement that one should consider 
an arbitrary puddle with N sites, and integrate over all N phases. Since the number of bonds exceeds the number of 
sites in two and three dimensions, it would be incorrect to consider the phase differences between bonds as independent 
variables. The maximum number of independent phase differences in a puddle of N sites is iV — 1. Starting from an 
arbitrary site in the puddle, a non self-intersecting string of — 1 bonds spans all N sites. The string, however, is 
not unique, hence in the partition function Zn for TV coupled sites, one must consider all possible self-avoiding string 
configurations of iV — 1 links. From the theory of self- avoiding random walks^^, it is known that the degeneracy Af of 
such configurations is 



{N-iy/^S^-\ d^2 

where Zd is an effective coordination number that depends on the dimensionality and the arrangement of grains. For 
a simple cubic lattice in three dimensions, Z3 = 4.68, slightly less than 6, which is the actual coordination number. 
Thus the contribution of an A-site puddle to the partition function, say in three dimensions, is 



(20) 




The optimum size of the puddle is reached when N — k, \/ in(7r/2g|i*7) ' ^^'^ dominant contribution to the 
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This result is valid under the condition 1 7 ^ 1/ y/f3Ec. 

We now have the necessary ingredients for calculating the conductivity a and tunnelling density of states Vi from 
Eq.(Pl and Eq.Q. Calculation of the conductivity a using Eq.© requires evaluation of a two-point phase correlation 
function Hj i+a, 

ftia+a = (cxp(-i(0i,i+a(T)-0i,i+a(O)))^. (22) 

The two points i and i + a should be chosen to lie in different puddles, for if they lie within the same puddle, 11 
would simply describe fluctuation of charge distribution inside a puddle; this contributes little to the conductivity a. 
This simplifies evaluation of the two-point phase correlation function to a product of two one-point phase correlation 
functions, ni,i+a ~ (exp(— i(^i(T) — (^i(0))))(exp(«(^i+a(''') — 0i-(-a(O)))). The averaging in Ea. l(^ should be performed 
over winding numbers {fci} as well as the phase fluctuations {(/>i}. The AES action in Eg. ljl^l) after integrating over 
the relative residual phase fluctuations then takes the form 



|i-j|=a 

+ E ^i^' cpiMM-u;,,) + ■■■ ■ (23) 



Upon performing the average, we obtain an expansion in increasing puddle size: 
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To calculate the conductivity a given by Eq.Q, we make the analytic continuation 17„ —iuj + e, and deform^'^ the 
contour of integration in the following manner: (0, (3) — > (0, ioo) + [ioo, ioo + P) + (ioo + (3, (3). For d.c. conductivity, 
we expand Eq.Q for small lo, and take the limit — )■ 0. Performing the integration yields the conductivity 
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Most of the contribution to Ea. (|25|l comes from two single-charge configurations (qNi,qN;+^) — (1,0), or (0,-1). 
In the former configuration, conductivity is dominated by (iVijA^i+a) = (1,A*), while in the latter configuration, 

conductivity is dominated by (A'^i, Ai+a) = (A*, 1), and N = \J \n(^^2k ^ usual. The result is 

(T - ^a^-"^ exp(-2 V/3£^c ln(7r/27id)). (26) 
We can similarly obtain the tunnelling density of states: 

^{^) ~ ^ cosh(/3£) exp(-2v//?£;cln(^/275d)). (27) 



III. CONCLUSION 



We propose that our simple model of a regular array may explain soft activation behaviour observed in real granular 
metala^'S'i'^. In real granular metals, inter-grain tunnelling may vary strongly between grains, but even in the presence 
of disorder, our physical mechanism could be applicable. Firstly, for weak disorder, suppose the inter-grain coupling 
for the i*'' tunnelling link has a distribution 7^ = 7-'^+'^' . Then for an A'^-site puddle, since the are random, the 



7 



tunnelling term Yiili 7^+Si«i ^ jg j^^^ seriously modified, and our conclusions hold. Secondly, as discussed 
in the context of granular superconductorsi^, theoretical calculations based on regular Josephson arrays seem to be 
relevant. The reason is that even for a wide distribution of couplings, only a narrow range of couplings is relevant, 
since (a) the extremely weak links can effectively be disregarded and (b) for links that are much stronger than 
average, one can approximate the connected grains as one single grain. While the tunnelling probability changes 
exponentially with length, the charging energy changes only linearly, so the variation of charging energies is relatively 
small. The system then effectively consists of such renormalised 'grains' linked by tunnelling of similar magnitude. 
Thirdly, if conduction occurs through a few ID paths, our result, being dimensionality independent, still applies. The 
observationsSi^, according to our picture, are robust even upon application of strong magnetic fields (> 10 tesla) and 
are independent of dimensionalitjiS. Nevertheless, further work needs to be done to understand properly granular 
metals with strong variation in inter-grain couplings. 

The AES approach, we use, views conduction as a Fermi Golden- Rule type incoherent tunnelling process. The 
obvious difference between our picture of soft activation and the Efros-Shklovski:^^(ES) theory, which also gives a 
similar temperature dependence of conductivity, is the on-site charging energy cost Ec in our model and lack of 
thereof in ES theory. Furthermore, the mutual interaction of charges (and excitonic effects) on widely separated 
grains plays no significant role in our analysis unlike in ES theory. Since the soft activation mechanism involves only 
nearest-neighbour hopping in comparison with long-distance variable range hopping, the magnetoresistance of soft 
activation here is expected to be very weak, which is consistent with experiments*. Another possibility^ considered 
in the literature suggests that the observed soft activation could be an artifact of a special distribution of grain sizes. 
Such a hope is belied by observation^ of the same soft activation in samples with a very narrow distribution of grain 
sizes. Also if we accept the conduction process as proceeding through tunnelling of charge between neighbouring 
grains, there would be little likelihood of finding the percolation paths in the wide range of temperatures through 
appropriately sized grains, should they exist, as neighbours. 

The relevance of our results as well as Refsi^*i^ should be explored beyond carefully prepared granular arrays. 
Recently, a logarithmic temperature dependence of conductivity in strong magnetic fieldsSSiSi and granular (or domain) 
structure^^ has been observed in certain underdoped cuprates. The insulating phase (even more underdoped) in the 
same materials exhibited the soft activation behaviour^'^ , which may be due to the mechanism proposed in this paper. 

In conclusion, our analysis of transport in granular arrays at not too low temperatures T ^ max {d, gd) in the 
framework of the AES approach shows that the transitions from a logarithmic temperature dependence of conductivity 
for strong inter-grain coupling (g 3> 1) to the soft activation behaviour for intermediate coupling (g ^ 1, 1 ^ 7 3> 
l/y/l3Ec) and further to the hard activation behaviour for weak coupling g <C 1 can be understood as arising from 
the competition between Coulomb blockade and tunnelling. This analysis is strictly valid for regular arrays and may 
be considered for experimental system s -^i -M . '? . 
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Appendix 

Note added in proof: The main idea of our paper is that charge sharing among several grains occurs at certain 
low temperatures, and the probability of sharing a charge between two grains is not an exponentially small function 
of temperature. In this Appendix we discuss various arguments (in addition to the calculation in the main text) in 
support of this main idea, which we developed after the original manuscript was submitted. We hope to publish a 
more detailed discussion elsewhere. 

We discuss in detail the situation of two connected grains to demonstrate again the essential physics. For two 
grains, the AES action can always be expressed in terms of the average phase (f)av — (f/*! + 02)/2 and relative phase 
(j) = ((/)i — 02), where (jfi and (/)2 are the phases on first and second grains. 
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The part of the action for the average phase is trivial and is easily transformed to the charge representation. In doing 
so it is necessary to satisfy carefully correct Matsubara boundary conditions of the original fields (pi and (l>2 ■ 

The crucial question is what the minimum charging energy of two grains is. Is the minimum charging energy Ec/2 or 
still Ec as for a single grain? For a two-grain system, this issue can be addressed by considering the phase correlation 
function for one of the grains, 

Ci(t) = {cos{Mr)-Mm 

= (cos((/)a.(r)-0a.(O)))(cos(((/)(T)-0(O))/2)). (29) 

The part of the action for the average phase 4>av corresponds to the charging energy Ec/2 for the total charge (gi +(72), 
quantised and equal to one. Prom this part, Ci(r) gets a contribution exp{—EcT/2). If g were zero, then the relative 
phase contribution to Ci(t) is also exp(— £'ct/2), so that Ci(r) = exp(— Sct), in accordance with our expectation 
for isolated grains. For a finite value of g, we demonstrate below that correlation functions of the relative phase 
fluctuations at large r decrease only algebraically as a function of temperature, and do not show a hard Coulomb 
gap2^. The gaplessness of the relative phase fluctuations unambiguously proves that the minimum charging energy of 
the two grains is halved, Ec/2, and associated only with the average phase. If instead, the relative phase correlator 
were gapped, with some effective charging energy E*, then the correlation function would decrease exponentially 
at long-r with the corresponding charging energy. Therefore the question about the effective charging energy is 
equivalent to the question of considering the long-r asymptotics (or equivalently low temperatures) of the correlation 
function Ci(t). The last statement is of a general character, since the long-time (or low temperature) asymptotics 
always reveals the lowest energy excitations (or configurations of the (/)-field) of the system. Namely, if the charge gap 
exists, this will become evident as an exponential decay of the corresponding correlation function at long times. In 
the literature, the charge gap is occasionally related to the amplitude of Coulomb blockade oscillations as a function 
of a gate voltage on a grain. In our case, the charge gap is the cost of putting one excess charge on a grain. These 
two definitions are not necessarily the same. For calculating the conductivity of the granular system, our definition 
of the charge gap is the appropriate one. 
Let us consider the relative phase correlator 

C{t) = (cos(0(r) - 0(0))), (30) 

which has been extensively studied in the literature, so a comparison with known results is possible. Besides, C(t) < 
(cos(((/)(t) — (f)(0)) /2), so if we can show that the large r behaviour of C(t) does not have a hard gap, then it is true 
for (cos((0(r) — (f)(0)) /2) too. We claim (and we are not the first ones) that the correlation function C(t) of the 
single-phase action decays, in fact, as a power-law, (T/T»)^/ sin^(7rTr), at very large r and not exponentially. Here T^, 
is an energy scale exponentially small in g. Thus the correlation function of relative phase fluctuations is not gapped 
for long times. This is a crucial point because a temperature dependence of C(t) that is not exponentially small 
in temperature at large r invariably leads to a soft activation behaviour of conductivity at low temperatures (see 
concluding remarks in this section) . Several arguments based on general results of statistical physics and mesoscopics 
as well as our numerical Monte Carlo simulations prove beyond doubt that the correlation function C(t) decays 
algebraically (proportional to 1/ sin^(7rTT)). First, note that the action for relative phase fluctuations is a one- 
dimensional field theory with a long-range interaction, gT^/ sin^(7rTT), in imaginary time. A general theorem of 
statistical physics due to Griffithg^^ states that the correlation function C(r) cannot decay faster than the interaction, 
gT^ / sin^ (ttTt) . The exponential decay is much faster than algebraic decay and therefore not possible. Second, it is 
widely recognised in the mesoscopics literature that at low temperatures the tunnelling to a quantum dot is dominated 
by so-called inelastic (or elastic) cotunnelling processes^SiSi. In the case of two grains, cotunnelling processes which 
are second order processes in the conductance g correspond to creation of electron-hole pairs on both grains. These 
processes are the lowest-energy gapless processes which can be closely associated with long-r behaviour of C(t), 
once again this demonstrates that the charge gap at low temperatures is effectively zero. In fact, the picture of the 
charge sharing can be equally well discussed in terms of the balance between Coulomb blockade and cotunnelling 
processes for a finite set of grains. Our results concerning cotunnelling processes are somewhat non-trivial, because 
we describe these processes in terms of the parameter j(T) for g 3> 1 (unlike the originally considered case oi g <^ 1 
of RefiSSiSi). Third, we undertook numerical simulations of the single-phase AES action using the path integral 
Monte Carlo method. It is possible to calculate directly by this method, without any approximations, not only the 
correlation function C(r) but also the parameter 7 as a function of g and T. Numerical results show clearly that 
the correlation function behaves as (T/T*)^/ sin^(7r Tr) in the large-r limit. Note that a rough estimate of C(t) as 
C(r) ~ EfcCOs(27rfcTT)exp(-5(fc))]/exp(-S'(fc)) using Eq.@ gives C(r) - sinh^ (7r7/2)/[sinh2 (717/2) -h sin2(7rrT)]. 
For small values of 7, the gaussian approximation is inaccurate, nevertheless the Lorentzian long time behaviour 
of C(t) inferred from Eq.Q clearly anticipates the (T/T*)^/ sin^(7rrT) result of exact numerical calculations, with 
7 ~ T/T* < 1 and large r. 
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In what follows we summarise the results of the path integral Monte Carlo simulations (the description of the 
method and further results will be published elsewhere?*). In FigQwe present the correlation function C(r) for 5 = 1 
and g — 1.5. One can see clearly that at short imaginary times-r, C(r) = 1 — l/{TTg) \n{gEc/2T), which is consistent 
with gren Calculated in Eq.Q ignoring winding numbers and various others'^. At long times-r, in Fig^ the correlation 
function decays as 1/ sin^(7rTT). 




FIG. 1; InC(r) vs Inr for g = 1 (red) and g = 1.5 (blue), at /9 = 768. The data (solid curves) show a crossover from logarithmic 
behaviour C"^"^^'"^ (dotted curves) to the power-law behaviour C^'^ (dashed). Discrepancies at small r are due to discretisation 
error. Note that C^^iog jg jjgt exactly a power law (or a straight line on a log — log plot): the cosec squared flattens out at 
r ~ (3/2. 



The detailed behaviour of the parameter 7 is given in Fig|21 as a function of g and in Fig|3| as a function of 
temperature T. The exact calculation of the determinant of the residual fluctuations in the gaussian approximation 
gives instead of 7r7|fci2| in Eq|51the following expression: 

7r7|/c| = Trg\k\ — 



F(l + fc + F (1 + fc + ^+^^^+^ 



F(l + A;)2 F(l + x) ' ^^^^ 

where x = and for brevity we denote the relative winding number fci2 = k as k. This expression is more precise 
than the expression (see Ea llO() calculated using Stirling's approximation. 
FigEl compares the derivative of the action, 

A^,{g) = d[S{k = 1,5) - S{k = 0,g)]/dg, (32) 

evaluated numerically, with various analytic approximations. The blue dotted curve is the gaussian approximation of 
Eq. and the blue dashed curve is Eq. (|1(J|I which can be shown to follow from Eq. using Stirling's approximation 
for the gamma functions. The red dashed curve corresponds to setting 7(T) = gren- 

Fig 13 compares 7(T), calculated numerically, with the gaussian approximation in Eq.lPU (blue dashed curve), and 
gren(T) Calculated from C(t), t oc /3, (top red curve) and gren{T) with an adjustable cutoff (lower red curve), as a 
function of temperature. 

We thus observe that the renormalisations of the quantities gren and 7 as a function of temperature are different 
for large g 1. Although this observation is not essential by itself for the charge sharing mechanism, the numerical 
simulation shows directly that the renormalisation of 7 is stronger (see below). The stronger renormalisation of 
7 in comparison with gren 

should presumably be associated with nearly zero-modes of residual fluctuations which 
exist around winding number trajectories. Namely, the gaussian fluctuations are stronger around winding number 
trajectories than around non- winding (fc = 0) trajectory because the square averaged fluctuations for zero modes (for 
n < k) are much stronger < 0^ >zm= 4:Ec/{ttT) than < cj)^ >~ 1/g for simple gaussian residual fluctuations. Note 
that zero-mode fluctuations need to be considered beyond gaussian approximation for Ec » T, because < >zm 
becomes easily much larger than (27r)^. Therefore fluctuations beyond logarithmic renormalisations are naturally 
expected and do occur as seen numerically. 

We end with two remarks. Since the temperature dependence of the relative phase correlator at large t is not 
an exponentially small function of temperature, but only a power law, optimising the probability of a charge shared 
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FIG. 2: Ail)i(g), defined as the derivative with respect to g of the normalised action (see Eg. 13211 1. plotted against g for 
/3 = 16. Solid black line: Monte Carlo. Dashed black line: small-g perturbation theory. Dotted blue line: Ai/jf^"''^, using 
the gaussian approximation of Ea. H^lfl . Dashed blue line: Aipf^^^^'^ , using Ea. ljlO^ which can be deduced from Ea. H^ljl using 
Stirling's approximation. Dashed red line: (1 — 1/ng), i.e., the derivative of grm with respect to g. 
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FIG. 3: 7 versus /3 for g — 32. Solid black line: Monte Carlo calculation. Dashed blue line: 7(r) in Eg. 1)31^ obtained by 
considering gaussian fluctuations. Top red line: gren{T) extracted from C(r) with r proportional to /?. Bottom red line: 
greniT) with adjustcd cutoff, gren{T) ~ [g — (I/tt) ln(const.i?c/T')]. Clearly, (l/7r)lnT is not in good agreement with Monte 
Carlo results for any adjustment of the cutofi'. The gaussian approximation gives a much better agreement with Monte Carlo 
results. 



among N grains, Pn ^ -f{T)^^^ exp{—Ec/NT) w exp{~2^/E^lIl(J-^^^)/T), gives a temperature dependence of the 
exponent that is always weaker than Arrhcnius' law, P ~ exp(— £'*/T). Thus the temperature dependence of the 
optimum probabihty P/v is the soft activation behaviour. Second, it has not escaped our attention that even for .g < 1, 
inelastic cotunnelling processes should make the charge sharing possible2&. 
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